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I. INTRODUCTION 

The geometrization of physical theories is a successful and challenging area in theoretical 
physics. The most well known examples are Hamiltonian mechanics based on symplectic 
geometry, General Relativity based on semi-Riemannian geometry and classical Yang-Mills 
theory which uses fibre bundles^. 

Geometric ideas have also found a clear utility in non-relativistic quantum mechanics 
problems because quantum theory can be formulated in the language of Hamiltonian phase- 
space dynamics^. Hence, the quantum theory has an intrinsic mathematical structure 
equivalent of Hamiltonian phase-space dynamics. However the underlying phase-space is 
not the same space of classical mechanics, but the space of quantum mechanics itself, i.e., 
the space of pure states or the space of mixed states. 

As opposite to what happens in General Relativity or in Gauge Theory where the metric 
tensor or the connection are related with the physical interaction, the most usual formulation 
of the geometry of non-relativistic quantum mechanics is not dynamic, in the sense that is 
insensitive to changes in the Hamiltonian of the system. Under these assumptions, that 
approach only makes use of the differential structure of the Hilbert space for quantum 
states and the Fubini-Study metric. See for example the geometric interpretation of Berry's 
phased. 

From a more dynamical point of view, A. Kryukov^ has stated that the Schrodinger 
equation^ for a pure state \at) (Plank's constant is set equal to 1) 



can be considered as a geodesic flow in a certain Riemannian manifold with an accurate 
metric which depends on the Hamiltonian of the system. 

The goal of this paper is to generalize the work of Kryukov for mixed states. To this 
end, we provide an underlying differential manifold to describe mixed states and a dynamic- 
dependent Riemannian metric tensor to analyze their temporal evolution. 

The mixed states are characterized by density matrices and the equation which plays the 
role of the Schrodinger one is the von Neumann equation^^ 




iH \at) 



(1) 



dpt 
dt 



(2) 
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To obtain the underlying differential manifold, following the Uhlmann's geometrization 
for non-relativistic quantum mechanics^Si"— , we make use of a principal fibre bundle such 
that the base manifold of that principal bundle is the space of mixed states. Finally, to 
provide the Riemannian metric we choose an appropriate metric in the principal bundle, in 
such a way that the projection of the geodesic flow in the principal manifold to the base 
manifold is just the temporal evolution given by the von Neumann equation. 

Among the geometric properties that are observed due to the movement of this geodesic 
flow, we analyze in this paper the phase volume conservation according to the Liouville 
Theorem. That allow us to show a geometric proof of the Poincare recurrence theorem 
relating it with the recurrence principle for physical systems with discrete energy levels. 
Let us emphasize that our geometric proof for the quantum Poincare recurrence is closer 
to the proof for classical mechanics^ (that also uses the conservation of the volume in the 
phase-space evolution) than the previous given in the quantum setting^i^i^. 

II. DENSITY MATRICES SPACE AS A BASE OF A PRINCIPAL FIBRE 
BUNDLE 

The most general state, the so-called mixed state, is represented by a density operator in 
the Hilbert space H. In this paper we assume always that dim('H) = n < oo, being Ti a 
vector space over the complex field ("H = C"). The density operator p is in fact a density 
matrix. Recall that a density matrix is a complex matrix p that satisfies the following 
properties: 

1. p is a hermitian matrix, i.e, the matrix coincides with its conjugate transpose matrix: 
P = P^- 

2. p is positive p > 0. It means that any eigenvalue of A is non-negative. 

3. p is normalized by the trace tr(p) = 1. 

Let us denote by V the space of mixed quantum states. Note that the space of pure 
states ViTi) is just 

Vin) = {peV\p' = p} . 

Recall that the space of quantum pure states has an elegant interpretation as a [/(l)-fibre 
bundle S{T-L) — )■ V{'H). Following Uhlmanrt^^"— and Bengson and Chruscihski books^^'^, we 
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can use a similar argument to that made in the case of quantum pure states. The key idea 
of Uhlmann's approach is to hft the system density operator p, acting on the Hilbert space 
7{, to an extended Hilbert space 

In quantum information theory-*--, the procedure of extension, T-i — t- is known as at- 
taching an ancilla living in "H. Obviously, the space of squared matrices A^„^n(C) (n rows, 
n columns) over C (that is a 2n^ real dimensional manifold) can be identified with 'H'^^^ 

As tr(PFTy'') is a smooth function over the space of squared matrices, by the Regular Level 
Set Theorem^, the set 

5o := {W e MnAQ ■ tT{WW^) = 1} , (3) 

is a smooth manifold of A^„ „(C). Actually, it is no hard to see that Sq is diffeomorphic 
to the sphere S^" -i. if 

p is a mixed state in V, we shall denote an element W G Sq a 

purification of p if 

p = WW^ , (4) 

therefore, we get the space of density matrices V by the projection n : Sq ^ V , where the 
projection is given by 

7r{W) = WW^ . (5) 
Observe that, if u is an unitary matrix (i.e, mm^ = m^m = I„) then 

7r{Wu) = n{W) . (6) 

Moreover, to fix notation recall that the Lie group U{n) is a Lie transformation group^^ act- 
ing on Sq on the right. In general, a principal fibre bundle^^ will be denoted by P{M, G, it), 
being P the total space, M the base space, G the structure group and vr the projection. For 
each X G M, n~^{x) is a closed submanifold of P, called the fibre over x. If p is a point of 
n^^i^x), then 7r"^(x) is the set of points {pa, a E G}, and is called fibre through p. 

At this point, an important question to answer, is if Sq{V ,U{n),Ti) is a principal fibre 
bundle over the base manifold V of density matrices. Unfortunately the answer is no, because 
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U{n) doesn't act freely on Sq. In general, Wu = W for W ^ Sq and u G U{n) doesn't imply 
that u = In, but observe that if det(Vr) 7^ ( i.e, W is an invertible matrix ) U{n) would 
act freely on our space. This is the way that we shall use to describe the space of density 
matrices. Instead of start with A^„^„(C), we start with the subset of invertible matrices. 
That subset has the differentiable structure of the Lie group GL(n,C). Then, we build a 
submanifold S of GL(n, C) given by 

S := {W eGL{n,C);tT{WW^) = 1} . (7) 

Finally, we obtain the base manifold using the projection vr : iS — ^ given by 

7r{W) = WW^ , (8) 

and therefore, S {V^ , U (n) , tt) becomes a principal fibre bundle. Observe that 

= {peV\p>0} , 

contains only strictly positive (or faithful) density operators. But V can be recovered from 
by continuity arguments^. 

In short, we describe the geometry of the density matrices as a base manifold of a principal 
fibre bundle consisting of a submanifold S of the Lie group GL{n, C) diffeomorphic to the 
sphere S^"^~^ as a total space and the Lie group f/(n) as structure group. 

Since S(V~^ , U{n), tt) admits a global section^° r : — )■ S 

rip):=VP , (9) 
therefore S{V~^ ,U{n),7r) is a trivial bundle from a topological point of view, that means 

S = V+ X U{n) . (10) 

III. HAMILTONIAN VECTOR FIELD, DYNAMIC RIEMANNIAN 
METRIC, SHG-QUANTUM FIBRE BUNDLE AND MAIN THEOREM 

In this section, we define a Riemannian metric for dynamics systems and we study how 
this metric acts within the tangent vector space of S. We also discuss the relationship 
that this metric has with other metrics such as the Bures metric or the metric proposed by 
Kryukov^. 
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A. Hamiltonian vector field, dynamic metric and its relation with other 
metrics 

In order to provide explicit expressions for tangent vectors to S and the metric tensor, 
we identify the tangent space TwMn,n{C) with A^„^„(C) itself. In such a way that since 
our total space 5 is a submanifold of the manifold A^„^„(C), where each point W & S is & 
matrix, and the tangent space Tw<S to S in the point W is a subspace of the tangent space 
TwAin,n{'C). We can use a matrix to describe a point W & S and a matrix to describe a 
tangent vector X G Tw<S too. 

First of all, note that the Hamiltonian operator H induces a vector field h : S ^ TS over 
S given by 

hw := -iHW , (11) 

where hw denotes the vector field in the point W E S, i.e, hw = h{W). That vector field h 
will be denoted as the Hamiltonian vector field. 

For any point W E S, and any two tangent vectors X,Y E TwS, we define the dynamic 
Riemannian metric gH{X,Y) as 

gniX, Y) := ^ti{X^H-'Y + Y^R-^X) . (12) 

It will be denote by V'^ the Levi-Civita connexion (the sole metric torsion free connexion) 
given by gn- In the definition f|T2l) we use assuming that H is an invertible matrix, 
but that in fact makes no restriction on the Hamiltonian of the system because we can set 
if — )■ if + I„ without change the underlying physics. Is not hard to see that g^ defines a 
positive definite inner product in each tangent space TwS, being therefore gu a Riemannian 
metric. 

With that metric tensor g^ the (sub)manifold {S,gH) becomes a Riemannian manifold. 
In order to fix the notation we denote {S{V^, U (n), tt), h, gn} the SHg-quantum fibre bundle 
of dimension n. 

The rest of this section will devote to study the inherited metric in the base manifold 
(theorem 1) from the the dynamic metric in the principal manifold and its relation between 
other metrics. 

The tangent space TwS at the point W E S can be decomposed in its horizontal Hw 
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and vertical Vw subspaces: 

Observe moreover that the vertical subspace Vw are the vectors tangent to the fibres. So, 
any vertical vector Xy G TwS can be written as 

Xv = WA , 

where A G u{n) (i.e, A is an antihermitian matrix). Note that our metric qh defines a 
natural connexion as follows: A tangent vector X at is horizontal if it is orthogonal to 
the fibre passing through W, i.e., if 

gHiX,Y)=0 , 

for all vertical vector Y at W {Y & Vw)- Hence X G TwS if horizontal if 

X^H-^W - W^H-^X = . (13) 
Therefore, we can define a metric gj^^ in the base manifold at any point p G V^, given by 

where Y,Ze TpV^ and Inor (respectively ^Hor) is the horizontal lift of Y (respectively Z). 
Theorem 1. The metric g^^ in the base manifold at any point p G can be obtained as 

gt{Y,Z):=Ur{H-'GyH-'Z) , 
where Gy is the unique hermitian matrix satisfying 

H-^YH-^ = GyH-^pH-^ + H-^pH-^Gy . 

Note that matrix Gy exists and is unique by the existence and uniqueness of the solution 
of the Sylvester equation^"^. Observe moreover that when H is the identity matrix 

gl\Y,Z):=Uv{GyZ) , 

where Gy is the (unique) solution of 

Y = Gyp + pGy , 

and that is the Bures metrio^^. 
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Proof. Let : R — )■ 5 be a curve, such that W is an horizontal vector, then 
Note, let us define A = H~^W so 

A^A = A^A . 

It is easy to see that the latter condition is fulfilled if 

A = GA , (14) 
where G is an Hermitian matrix. Therefore 

W = HGH-^W . 

Hence applying equation H] 

7r,{W) = WW^ + WW^ = HGH-^p + pH-^GH . (15) 

Suppose that 

TT,{W) = Y n,{V) = Z 
W{0)W{Oy = V{0)V{Oy = p , 

then 

gt{Y, Z) =gH{W, V) = hi{W^H'^V + V^R-'^W) 

= -ir{H-^GYGzH~^p + R-^GzGyH-^p) , 

where 

W = HGyH^^W V = HGzH-^V . 
Applying equation f|T5l) in n^iV) 

Z = HGzH-^p + pH-^GzH . 
Using the above expression (fT6!) the theorem follows. □ 

In the case of pure states, our Hilbert space is C" and the tangent space will be C"' 
too. Following Kryukov^, we can define a metric gK{X,Y) for any two tangent vectors 
X = {x,x*), Y = {y,y*), by 

gKiX,Y) ■.= Re{{H-'X,H-'Y)) , 

where {X,Y) = Yl^=i^iyh therefore 

gK{X, Y) := i {{H-'X, R-^Y) + {R-^Y, R-^X)) = hi^X^R^^Y + Y^R-^X) . 

When R is the identity, we recover the Fubini-Study metric. 
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B. Geometric structure of the SHg-quantum fibre bundle 

As we have seen in the previous sections in the SHg-quantum fibre bundle {S{V'^, U (n), tt), h, gn} 
of dimension n, S{V^ iU{n),'K) is a principal (and trivial) fibre bundle, S is diffeomophic 
to the sphere of dimension — 1, h is & vector field on 5, and {S^qh) is a Riemannian 
manifold. But the SHg-quantum fibre bundle has more geometric properties : 

Theorem 2 (Main theorem). Let {S{V^, U{n),T^), h, qh} he a SHg-quantum fibre bundle of 
dimension n. Then: 

1. h is a Killing vector field of [S^gn)- 

2. The integral curves 7:/c]R— >-<S of h are geodesies of {S, gn) ■ 

3. The projection on the base manifold of the geodesic 7 satisfies the von Neumann 
equation 

— TT O 7 = — i [if, TT O 7] . (17) 

at 

Proof. Condition (1): In order to proof that /i is a Killing vector field, we only have to show 
that the fiow : <S — >■ 5 given by 

^JW) = PV, where W eS 

(18) 

is an isometry, i.e, for any X,Y & T^S 

gH{vu{X):Vu{Y))^gH{X,Y) . (19) 

Note that 

^t,(X)^e-'"'X , (20) 

and 

gHi^UX), ^UY)) =gH{e-'"'X, e^^V) 

=-tr {{e-''''X)^H-\-'"'Y + {e-'"'Y)^ H-^e-'^^'X) 

I (21) 

= -tr (Xte^^*i/-2g~.my ^ Y^e'^'H-^e-^'^'X) 

=^tr (Xt//- V + Y'^H-'^X) = gniX, Y) . 
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Conditions (2) and (3): First of all observe that if 7 is the integral curve of the vector 
field /i, i.e, 

^ = h^ = -iH-f . (22) 

The projection of 7 satisfies 



j^Alit)) = I (7(^)7^^)) = 7(^)7^^) + 7(^)7^^) = 7(^)7^^) + lit)my 

(23) 

= -iH^^\t) + -f{t){-iH-fy = -i[H, 7r(7(t))] . 

Hence, the projection of the integral curves of the vector field h satisfies the von Neumann 
equation. So all we have to prove is that curves are actually geodesic curves 



Vf,/^, = . (24) 

Since h is a. Killing vector field, we only have to proof that h is an unitary vector field (due 
any unitary Killing vector field is a geodesic). Namely, the equality 

gnih,, h,) =gHHHj, -iHi) = tr {{-iH^Y R-^-iH^)) 

(25) 

=tr (j^HH-^Hj) = tr (7^7) = 1 . 
Finally, since "H is an unitary Killing vector field and the integral curves of any Killing vector 
field of constant length is a geodesic (see appendix theorem |H]), the integral curves of "H are 
geodesies. □ 

Remark. Let us emphasize that for any hermitian matrix A = , we can build the vector 
field Aon S given by —lAW for any W E S. It is easy to check in a similar way to equation 
[21] that if [if, A] = 0, ^ is a Killing vector field. Therefore, the set of operators compatibles 
with the Hamiltonian are related to the set of isometrics of (iS, qh)-, and we can identify any 
conserved quantum observable with a Killing vector field. 



IV. GEOMETRIC APPROACH TO QUANTUM POINCARE 
RECURRENCE 

As we know from the main theorem, /i is a Killing vector field on the principal manifold 
{SjQh) endowed with the dynamic metric gn- Then, the transformations given by the 
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1— parametric subgroup (pt '■ S ^ S of integral curves of h are distance-preserving and 
volume-preserving (see appendix theorem [9]). That two facts directly have the following 
consequences 

Theorem 3 (Insensitivity to Initial Conditions Theorem). Let {S{V'^ ,U{n),7T), h, qh} be a 
SHg-quantum principal bundle of dimension n. Then, for any two points W,V E S 



being the (ft the 1— parametric subgrup of transformations given by the integral curves of the 
Killing field h. 



Figure 1. Since the Hamiltonian vector field is a Killing vector field, its flow (p preserves the volume 
(Lioville theorem) in the sphere 5, where each point can be projected into the space of density 
matrices . 

The classical Liouville theorem^' states that the natural volume form on a symplectic 
manifold is invariant under the Hamiltonian flows. In our case, we have the 1— parametric 
subgrup of transformations (ft'-<S^S given by the integral curves of the Killing fleld h 
and we can set 



dist{iptiW),iptiV)) = dist{W,V) 



(26) 




GL(n,C) 
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Theorem 4 (Liouville Type Theorem). Let {V~^ ,U{n),'n'), h, qh} be a SHg-quantum prin- 
cipal bundle of dimension n. Then for any domain Q G S 

Vo\{ipt{^)) = Vol(fi) , (27) 

being the ipt the 1— parametric subgrup of transformations given by the integral curves of the 
Killing field h. 

Using the above theorem, we can therefore state a similar theorem to the Poincare recur- 
rence theorem^. 

Theorem 5 (Poincare Type Theorem). Let {S{V^ ^ U (n), vr), /i, qh} be a SHg-quantum prin- 
cipal bundle of dimension n. For any domain Vt d S and any time period T G there 
exist a point x and a positive integer k > such that 

cpkTix) e n , (28) 

being (ft'-S^S the 1— parametric subgrup of transformations given by the integral curves 
of the Killing field h. 

Proof. Consider the following sequence of domains 

VL, ipri^), v?2r(fi), ■ ■ ■ , Vkri^), ■ ■ ■ 

All domain in the sequence belongs to the same volume Vol(fi). If the above domains never 
intersect S would obtain an infinite volume, but S is compact, so Vol(5) < oo. Then, there 
exist / > and m > I such that 

^lT{n)nipmT{^) , (29) 

so 

nn^^^.i)T{^)j^iD . (30) 

Setting k = m — I the theorem is proven. □ 
Joining the above theorem with the Insensitivity to the Initial Conditions we get 

Theorem 6 (Strong Poincare Type Theorem). Let {S{V^ ,U{n),TT), h, qh} be a SHg- 
quantum principal bundle of dimension n. Then, for any point W & S, any e > and any 
T G M^, there exist a positive integer k > such that 

dist{W,ipkT{W))<e , (31) 
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being (ff-.S^S the 1— parametric suhgrup of transformations given by the integral curves 
of the Killing field h. 

Proof. Let us consider the domain 

B^{W) = |v e 5 : dist {W, y) ■ (32) 

Applying now the Poincare type theorem there must exist Wq G B^{W) and k > such 
that 

cpkTiWo)eB,^{W) . (33) 

So, 

dist {W,ip,T{Wo))<^ . (34) 
But, by the Insensitivity to Initial Conditions Theorem 

dist {ifiMW), ipkTiWo)) = dist {W, Wo)<^ . (35) 

Therefore, applying the triangular inequality 

dist (W, ifikAW)) < dist (W, ifikxiWo)) + dist (ifikAWo), ^kxiW)) < t . (36) 

□ 

A. Physical systems with discrete energy eigenvalues 

Using previously stated theorems we can give an alternative proof and more geometric 
sense of well-known^i^'^ principle of recurrence for physical systems with discrete energy 
eigenvalues. 

Thus, defining the length ||A|| of a matrix A as follows^ 

pl| = Vte(lM) . 

Then 

Theorem 7. Let p he a mixed state of a quantum system with discrete energy spectrum. 
Then, p is almost periodic. Namely, for an arbitrarily small positive error e the inequality 

\\p{t + T) - p{t) II < e for all t (37) 
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is satisfied by infinitely many values of T, these values being spread over the whole range 
—oo to oo so as not to leave arbitrarily long empty intervals. 

Proof. Let p{t) be the density matrix of a system with a discrete set of stationary states, la- 
beled n — 0,1,2, ■ • ■ , with energies En, some of which may be equal if there are degeneracies. 
In energy representation the matrix elements are 

Pnn'it) = {n\p(t)\n') . (38) 

Let Tn — \n){n\ be the projection operator onto the nth stationary state, then 

p'^"'(t) = r„p(t)r„, , (39) 

is the matrix which in energy representation has only one nonzero element, equal to Pnn'{t) 
and in the location {n,n'). These matrices are orthogonal in density space 

(p""'(i),p-"""'(t)) = 6nn"6n'n"'\pnn'{t)\^ , (40) 

and 

oo oo 



n=0 n'=0 

oo oo 

n=0 n'=0 

where ujnn' = (-E'n' ~ En). Now, consider the finite sum 

JV N' 



(41) 



a--'it) = ^^r'{t) (42) 



n=0 n'=0 

as an approximation to p{t). The square of the error is 



Mt)-a^^'{t)r = \\ y: e p'"' 

n=N+l n'=N'+l 



E E wp'^'^'itn" (43) 

=Af+l n'=N'+l 
oo oo 

E E ii^'^"'(o)ii' 



n=N+l n'=N'+l 



The second equahty follows from the orthogonahty of the p""' . Since the error is indepen- 
dent of the time, a^^' (t) converges uniformly to p{t) (in the || ||-norm sense). So p{t) can 



14 



Geometric approach to Quantum Dynamics 



be approximate by a^^'{t). Since a^^' (t) is a discrete density with finite energy levels, 
a^^'it) e P+, and the set 

(^^^') := {p G P+ : Hp - a^^'m < e} , (44) 

is an open precompact set in . Using the global section given in equation (|9]), t{B^) will 
be an open precompact set of S. But applying the Strong Poincare Type Theorem for any 
time period T > there exists k > such that 

dist(r(a^^\t)),^,T(r(a^^'(t)))) =dist(r(a^^'(t)),r(a^^'(t + A:T))) <e , (45) 

for any e > 0. Namely, 

r(a^^'(t + fcT))Gi?f(r(a^^'(t))) , (46) 

being Bf{T{a^^'{t))) the geodesic ball in S centered at T{a^^' (t)) of radius e. Now choosing 
e small enough 

i?f(r(a^^'(t)))Cr(i?r(a^^'(t))) . (47) 

Therefore by (l46ll 

r(a^^'(t + kT)) e r(5r(cr^^'(t))) • (48) 
Projecting to the base manifold 

cr^^'(t + fcT)e5r(cT^^'(t)) . (49) 
And, by definition of Bf^{a^^'{t)) 

||(T^^'(t + A;T)-a^^'(t)|| <e . (50) 

And the theorem is proven. □ 

V. APPENDIX 

In this section we recall several well known results about Killing fields on Riemannian 
manifolds, for a more detailed approximation see^i 

Theorem 8. (see also^) Any integral curve 7:JcM— i-M of a Killing vector field X of 
constant length sJg^X, X) is a geodesic on {M,g). 
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Proof. Here, we need 

V^7 = VxX = (51) 

but, since X is a Killing vector field the Lie derivative of the metric is zero Lxg = and 
(see^^ proposition 25) VX is skew-adjoint relative to g, so 

g{VxX, W) + g{VwX, X) = (52) 

for any X G TS. So 

= g{VxX, W) + l/2W{g{X, X)) = g{VxX, W) (53) 

then VxX = 0. □ 

Theorem 9. Let (M, g) be a Riemannian manifold, let X he a Killing vector field over M , 
and denote by (ft '■ M ^ M the 1— parametric subgroup of transformations given by X (i.e, 
Vo{p) = P> i^t{p)\t=o = Xp), then 

1. Given any two points p,q E M, dist{p,q) = dist{{pt{p) , ipt{q)) , 

2. Given any domain Q C M, 'Vol{(pt{^)) = Vol(r2). 

Proof. Let (fti^) be the flow of the domain Q. Allow us denote 

V{t):=Vol{ipt{n)) . (54) 

Then, the divergence is just (see^) 

V'{0) = [ div ndfXg^ , (55) 



n 



where dfig^ denotes the Riemannian density measure. 
The divergence is define as^ 

div H = tr(F ^ VyX) . (56) 
Given an orthonormal base in T^S 

v^n = J2 Y'K'H = E "^'anivln, e,)e, , (57) 
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where := gH{Y,Ei). Therefore 

diY n = J2gH{Vln,Ej) , (58) 

i 

But since H is a Kilhng field V^H is skew-adjoint relative to qh (see^^ proposition 25 
again), so 

div 7/ = . (59) 
And the theorem is proven. □ 
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